Astronomy & Astrophysics manuscript no. 6673 


© ESO 2008 


February 5, 2008 





o 
o 

(N 



in 



(N 

> 
o 

(N 
(N 



O 
O 



Morphological evolution of z ^ 1 galaxies from deep K-band AO 
imaging in the COSMOS deep field * 

M. Huertas-Company^'^, D. Rouan^, G. Soucail^, O. Le Fevre^, L. Tasca^, T. Contini^ 

' LESIA-Paris Observatory, 5 Place Jules Janssen, 92195 Meudon, France 

e-mail: marc .huertas@obspm. fr 
' Laboratoire d'Astrophysique de Toulouse-Tarbes, CNRS-UMR 5572 and Universite Paul Sabatier Toulouse III, 14 Avenue Belin, 

31400 Toulouse, France 

' LAM-Marseille Observatory, Traverse du Siphon-Les trois Lues BP8-13376 Marseille Cedex 12, France 
lAA-C/ Camino Bajo de Huetor, 50 - 18008 Granada, Spain 

Received jdate^ / Accepted idate,., 

ABSTRACT 

Context. We present the results of an imaging program of distant galaxies (z ~ 0.8) at high spatial resolution (~ 0.1") aiming at 
studying their morphological evolution. We observed 7 fields of 1' x 1' with the NACO Adaptive Optics system (VLT) in Ks (2. 16/jm) 
band with typical V ~ 14 guide stars and 3h integration time per field. Observed fields are selected within the COSMOS survey area, 
in which multi-wavelength photometric and spectroscopic observations are ongoing. High angular-resolution K-band data have the 
advantage of probing old stellar populations in the rest-frame, enabling a determination of galaxy morphological types unaffected by 
recent star formation, which are more closely linked to the underlying mass than classical optical morphology studies (HST). Adaptive 
optics on ground based telescopes is the only method today for obtaining such a high resolution in the K-band, but it suffers from 
limitations since only small fields are observable and long integration times are necessary. 

Aims. In this paper we show that reliable results can be obtained and establish a first basis for larger observing programs. 
Methods. We analyze the morphologies by means of B/D (bulge/disk) decomposition with GIM2D and C-A (concentration- 
asymmetry) estimators for 79 galaxies with magnitudes between = 17-23 and classify them into three main morphological 
types (late type, early type and irregulars). Automated and objective classification allows precise error estimation. Simulations and 
comparisons with seeing-limited (CFHT/Megacam) and space (HST/ACS) data are carried out to evaluate the accuracy of adaptive 
optics-based observations for morphological purposes. 

Results. We obtain the first estimate of the distribution of galaxy types at redshift z ~ I as measured from the near infrared at high 
spatial resolution. We show that galaxy parameters (disk scale length, bulge effective radius, and bulge fraction) can be estimated with 
a random error lower than 20% for the bulge fraction up to = 19 (AB = 21) and that classification into the three main morpho- 
logical types can be done up to K, = 20 (AB = 22) with at least 70% of the identifications correct. We used the known photometric 
redshifts to obtain a redshift distribution over 2 redshift bins (z < 0.8, 0.8 < z < 1.5) for each morphological type. 

Key words, galaxies: fundamental parameters - galaxies: high redshift - techniques: high angular resolution 

1. Introduction in particular the process and history of mass assembly. One 

, , , • , • classical observational way to test those models is to classify 

The process of galaxy formation aiid the way galaxies evolve is j^^^^^ accordi ng to morphological criteria define d in the 

still one of the key unresolved problems in modern astrophysics. ^^^^ ^j^^^^^^ (lHubblelll936t Ide VaucouleurslfT948t Is^i^ 

In the currendy accepted hierarchical picture of structure ^^.^^ ^^1^^^^ physi cal properties, and to 

formation, galaxies are thought to be embedded in massive dark foUpw this classi fication across time. (lAbraham et al.l [T996t 

^ISith^tawfrmn^nsity fluctuations in the early universe Lsjmard et aLlMlAbraham et al.ll20()3h . Indeed, since galaxies 

(| Fall&Efstafliiou| | l980|) and initially contain baryons in a hot ^^^^ recognized as distinct physical systems, one of the main 

gaseous phase. This gas subsequentlicools, and some fraction extragalactic astronomy has been to characterize 

exentual y condenses into stars^ ^ LiUy et al. | 1996 ; MadajietalJ their shapes. Comparison of distant populations with the ones 

[19981). However, many of the physical details remam uncertain, ^^^^^ ^^^^^ Universe might help to clarify the role 

* Based on ESO observations at the VLT Programmes P73.A-0814A "^^'^^^ °f "^^i" "^''"^'^ g^l^'^y evolution. 

andP75.A-0569Aandl75.A-0839 (zCOSMOS). (Cole et al. 2000; Baugh et al. 1996). Progress m this field 

** Based on observations obtained with MegaPrime/MegaCam, a joint ^as come from observing deeper and larger samples, but also 

project of the CFHT and CEA/DAPNIA, at the Canada-France-Hawaii from obtaining higher spatial resolution at a given flux and at 

Telescope (CFHT) which is operated by the National Research Council a given redshift. In the visible, progress has been simultaneous 

(NRC) of Canada, the Institut National des Science de I'Univers of the on those two fronts, thanks in particular to the ultra-deep HDF 

Centre National de la Recherche Scientifique (CNRS) of France, and fields observed with the Hubble Space Telescope. HST imaging 

the University of Hawaii. This work is based in part on data products has brought observational evidence that galaxy evolution is 

produced at TERAPIX and the Canadian Astronomy Data Centre as diff^erentiated with respect to morphological type and that a 

part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- j fraction of distant galaxies have peculiar morphologies 
rative project of the NRC and CNRS. ^ 6 f f 6 
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that do not fit into the elliptical-spira l Hubble sequence . 
(iBrinchmannet al.lfT998t IWolfet al.ll200l lllbertet al.ll2006bl) . 
These results can however be biased by the fact that most of the 
sampled galaxies are at large redshift and are analyzed through 
their UV rest-frame emission, which is more sensitive to star 
formation processes and to extinction. Moreover, it now seems 
clear that evolution strongly depends on the galaxy mass in 
the sense that massive systems appear to have star formation 
histories that peak at higher redshifts, whereas less massive 
systems have star formation histories that peak at progressively 
lower redshifts and ar e extended over a longer time interval 
(downsizing scenario: | Cowie et al. 1 119961; iBrinchmann & EUisI 
120001: lBundvetal.ll2005h . This~could explain the fact that the 
population of massive E/SO seems to be in place at z ~ 1 and 
evolve passi vely towards lower redshifts dPe Lucia et al. lHooi; 
IZucca et"al1 12006). However, most of these studies are based 
on spectral type classifications, and their interpretation in the 
framework of galaxy formation is not straightforward since 
galaxies move from one spectral class to another by passive 
evolution of the stellar populations. 

In this context, high-resolution near-infrared observations 
are particularly important because the K-band flux is less 
dependent on the recent history of star formation, which peaks 
in the UV in rest frame, and thus gives a galaxy type from 
the distribution of old stars that is more closely related to the 
underlying total mass than optical observations. This is why a 
large number of K-band surveys have been carried out using 
ground-based telescopes with different spatial coverages and 
limiting magnitudes in order to perfor m co smologica l tests 
by means of galaxy counts essentially dGardner et al.l 119931: 
Glazebrook et alj \l99^ iMcLeod et al.l 119951: iBershadv et alj 



19981: lMcCrackenetal.1120001: iMaiharaet alJ (20011) . However, 
no morphological information can be found in particular due to 
the seeing limited resolution, even with superb image quality as 
in the ongoing CFHT/WIRCAM survev. IConselice et alj (|2000|) 
prove that only when the ratio of the angular diameter substend- 
ing Q.SIi^^kpc at a given distance to the angular resolution of the 
image is around 1 can reliable morphological estimators such as 
the asymmetry be obtained. 

Consequently, adaptive optics (AO) installed on ground- 
based telescopes will be a powerful method for obtaining near- 
infrared high-resolution data in the future and an excellent com- 
plement to space observations, such as those that will be per- 
formed with the HST (WFC3, 2008) and with the JWST (2014). 

The use of classical AO for deep surveys suffers, however, 
from inherent limitations such as the non-stationary PSF on 
long integration times and the finite isoplanetic field. This is 
why preliminary studies to probe the accuracy of adaptive optics 
are required, before launching very wide surveys. In particular, 
it is important to determine whether automated morphological 
classifications can be performed. Indeed, given the large number 
of galaxies, such automated methods for morphological classifi- 
cations are desirable. 



iMinowa et al.l ( l2005h published first results based on 
adaptive optics observations. They achieved the impressive 
limiting magnitude of K ~ 24.7 with the Subaru Adaptive 
Optics system with a total integration time of 26.8/ir over one 
single field of 1' x 1'. They proved that the use of adaptive 
optics significantly improves detection of faint sources but did 
not ob tain morphological information. In a recently submitted 
paper, ICresci et"aLl (l2006h perform a morphological analysis 



based on AO data for the first time. They observed a 15 arcmiri^ 
in the band with NACO (SWAN survey) and classified 
distant galaxies into two morphological bins (late type, early 
type) by performing a model fitting with a Sersic law. They 
compared number counts and size-magnitude relation, for 
early and late-types separately, with hierarchical and pure 
luminosity evolution (PLE) models, respectively. They conclude 
that hierarchical models are not consistent with the observed 
number counts of elliptical galaxies and that PLE m odels are 
prefer red. However, as discussed in several studies (iGardnerl 
1 19981) . despite galaxy counts still being one of the classical 
cosmological tests, their interpretation remains diflicult. It is 
thus not realistic to expect galaxy counts alone to strongly 
constrain the cosmological geometry or even to constrain galaxy 
evolution. A more complete study needs redshift estimates, 
which is lacking in the SWAN sur vey. That is the main reason 
for h aving select the COSMOS dScoville & COSMOS TeamI 
I2005h field to conduct our pilot program since multi-wavelength 
photometric and spectroscopic observations are performed. This 
ensures a reliable redshift estimate for all our objects. 

In this paper, we continue this AO validation task by mor- 
phologically classifyi ng a sample of 79 ga laxies observed using 
param etric (GIM2D, 
(C-A, [Abraham et al 



Simard et al.l (l2002h ) and non-parametric 
1994 ) methods and comparing them. 



Fields, observed with NAOS/CONICA adaptive optics system, 
are distributed over a larcmin^ area. We obtain for the first time 
an estimate of the distribution of galaxies in three morphological 
types (E/SO,S,Irr) at redshift z ~ 1 as measured from the near 
infrared at high spatial resolution. We then use the photometric 
redshifts to look for evolution clues as a function of morphology. 

The paper proceeds as follows: the data set and the reduc- 
tion procedures are presented in the next section. In § [3] we fo- 
cus on the detection procedures and the sample completeness. 
In § |4] we discuss the estimate of redshifts using the multi- 
wavelength photometric data from COSMOS. The morphologi- 
cal analysis is discussed in §|5]using bulge/disk decompositions 
and concentration and asymmetry estimators. Simulations for er- 
ror characterization are carried out for both classification meth- 
ods, and comparisons between classifications are shown. In § |6] 
we compare the data with ground and space observations and 
use those comparisons to discus the morphological evolution in 
the last section. Throughout this paper magnitudes are given in 
the Vega system. We use the following cosmological parameters: 
//o = 70 km S-' Mpc-' and (Qm, ^k) = (0.3, 0.7, 0.0). 

2. The data set 

Seven fields of 1' x 1' were observed in Ks band (2.\6^m) 
with the NAOS/CONICA-assisted infrared camera installed on 
the VLT. Q The fields were selected within the COSMOS sur- 
vey are^ In order to ensure a reliable AO correction, relatively 
bright stars {V ~ 14) were selected. We added a color criterium 
(B - R < 1 .0) in order to benefit from a large attenuation of the 
flux in the near-IR and thus a smaller occultation of the central 
region of the images. The pixel scale (0.054") was chosen to be 
twice the Nyquist-Shannon requirements with respect to the tele- 
scope diffraction limit in order to have larger fields. With such a 
pixel sampling, the PSF FWHM would only be one pixel width 
in the limit of perfect AO correction. However, we will show in 



P73.A-0814A and P75.A-0569A 

,http : //www ■ astro . caltech . edu/cosmos/ | 
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Field 


a 


6 


exp.time (s) 


Seeing (arcsec) 


STARl 


10:00:16 


+02:16:22 


10350 


0.09 


STAR2 


10:00:52 


+02:19:52 


7650 




STAR3 


10:00:10 


+02:06:08 


7650 


0.12 


STAR4 


09:59:52 


+02:05:00 


7170 


0.08 


STARS 


10:00:14 


+02:09:09 


10200 


0.09 


STAR6 


10:00:02 


+02:06:57 


7650 


0.13 


STAR 10 


09:59:56 


+02:04:07 


10000 





Table 1. Summary of observations for the seven analyzed fields. 
The mean seeing is estimated when faint stars were detected. 




Radial Distance (arcsec) 

Fig. 1. Radial profile of all detected stellar objects with Ks < 
19. The mean resolution is ~ 0.1", broader than the telescope 
diffraction limit, as a consequence of the long integration time 
(~ 3/1) and the large fields. 



Sect. 5 that our PSF reconstruction procedure can handle such 
undersampled data. We also note that only partial AO correction 
was actually achieved during our observations, and our PSF re- 
construction is thus more than adequate here. Figure[T]compares 
the radial profiles of 5 detected stellar objects. Indeed, this pro- 
gram is pushed to its limits in terms of field size, exposure time, 
and brightness of the guide star 

Data are reduced in a standard way: exposures are taken in 
the auto-jitter mode, which means that the pointing is randomly 
shifted within a 7"box, in order to improve flat fielding, bad 
pixel correction, and sky background withdrawal. The sky value 
in each pixel is estimated by performing a clipped median of all 
the exposure frames: the 10% faintest and brightest values are 
removed from the computation. Cosmic ray and flat corrections 
were applied, and recentering was done using a standard cross- 
correlation method. Recentering is done at a sub-pixel level. For 
that purpose a cubic interpolation of resampling was performed. 
After stacking, a global estimate of the sky background was per- 
formed by computing the stacked image spatial median. The fi- 
nal image was obtained after substraction of this value. 




Fig. 2. Objects classification using flie MAG_AUTO-]V[UJVIAX 
plane. 

Photometric ze ropoints were determined using 2MASS stars 
einmannl[l992h . We performed aperture photometry on the 
guide stars and compared it to 2MASS data to deduce the zero- 
points. Note that the change of detector between periods P73 and 
P75 resulted in different zeropoint values for each run. The aver- 
age zeropoint for the first period is: 22.82 + 0.06 and for the sec- 
ond period 23.29 ± 0.06. W e also used the ESO calibration data 
set standard stars ( jPerss on et al. 1998) to validate our measure- 
ments (ESO pipeline computations and our own measurements 
on the standard stars images). There is good agreement between 
all these values. 



3. Detection and completeness 

All objects having a 3cr signal above sky, over 4 four contiguous 
pixels are detected using SExtractor (Bertin & Arnouts 1996|). 
We decided to apply this low detection threshold, even if the 
main goal of this paper is to perform a morphological analy- 
sis, for two main reasons. First we want to test the ability of 
adaptive optics based observations to obtain morphology, so we 
are seeking the limits; and second, we wanted to be sure that 
no objects are lost when computing number counts to compare 
with other near-infrared observations (see Appendix A). We find 
285 objects over the 7 fields. We then performed a cleaning 
task in order to separate galaxies from stars and spurious de- 
tections. This was made using the SExtractor MU_MAX and 
MAG_AUTO parameters that give the peak surface brightness 
above the background and the Kron-like elliptical aperture mag- 
nitude, respectively. The distribution of objects in this param- 
eter space clearly defines three regions that separate extended 
sources from point-like or non-resolved sources and from spuri- 
ous detections (Fig. |2]i. In this separation scheme, objects with 
very faint magnitude and high peak surface brightness are con- 
sidered as false detections. Boundaries were drawn manually and 
a visual inspection confirms that known stars in the field are in- 
deed identified as point sources. We consequently identified 79 
galaxies, 19 stars (or unresolved objects), and 187 spurious ob- 
jects in the whole field. 

The sample completeness for point sources was estimated by 
creating artificial point sources from fields stars (see sect. 5 for 
detailed explanations) with apparent magnitudes ranging from 
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Fig. 3. NAOS/CONICA K,-hand image of the field centered at a = 10:00:16, 6 = +02:16:22 . The total integration time is 10350 
s. The field size is 1' x 1' with a pixel scale of 54 mas. Circles are detected galaxies and boxes are stars. The stellar FWHM is 
measured to be 0.1". The bright star at the center of the image is used as the AO guide star. 



Ks - l?i to Ks - 24 and placing them at random positions. We 
ran SExtractor with the same configuration as for real sources 
and looked for the fraction of detected objects. We find that the 
sample is 50% complete at Ks - 22.5 (or AB = 24.5) for point 
sources. Completeness for extended sources is estimated in a 
similar way: we generate 1000 galaxies with exponential and 
de Vaucouleurs profiles of different morphological types (bulge 
fraction ranging from to 1) and with galaxy parameters uni- 
formly distributed. In particular, the sizes of disks and bulges are 
distributed uniformly between < rj < 0.7 and < < 0.7 
as detailed in Sect. § |5] This leads to half-luminosity radii rang- 
ing from < ri, < I . We find that the sample is 50% complete 
at Ks - 21.5 (or AB - 23.5) for this population of extended 
sources (Fig. |4]i. We used this completeness to compute num- 
ber counts and to compare it to other near-infrared surveys in 
appendix ??. 

4. Photometric redshifts 

Galaxy number counts provide a useful description of galaxy 
populations but suffer from numerous degeneracies when trying 
to trace the evolution of galaxy populations. Model predictions 
are subject to uncertainties in the spectral energy distributions 
and evolution of galaxies and in the free parameters specifying 
the luminosity function, the cosmological geometry, the number 
and distribution of galaxy types, and the effect of dust and merg- 
ing. The need to have redshift information is therefore the rea- 
son for selecting the NACO fields within the ongoing Cosmic 
Evolution Survey (COSMOS) in which multi-/l and spectro- 
scopic observations are performed. COSMOS is designed to 



17 . M 22 
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100% 50% 0% 

Fig. 4. Completeness for extended sources. Galaxies with pa- 
rameters {rj,r^,BIT) uniformly distributed are simulated and 
placed at random positions in the fields. 



probe the correlated evolution of galaxies, star formation, active 
galactic nuclei (AGN) and dark matter (DM) with large-scale 
structure (LSS) over the redshift range z = 0.5 to 3. The survey 
includes multi-wavelength imaging and spectroscopy from X- 
ray to radio wavelengths covering a 2 square deg area, including 
HST imaging of the entire field. 

All these data are used for a direct estimate of the photo- 
metric redshifts of the galaxies detected in the NACO fields. 
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using asymmetry estimators, and then we separated the regular 
objects between early and late type objects. 

Throughout this section, we use extensive simulations for 
error estimates and calibration of the automated classifications 
as explained below. For all the simulations, we assumed that 
bulges are pure de Vaucouleurs profiles (n - 4) and that disks 
are exponential profiles. We then generate galaxies with param- 
eters uniformly distributed in the following parameter space: 
< B/T < 1, 0" < rd < 0.7", O" < < 0.7", 0° < / < 70°, 
and < e < 0.7. Both bulge and disk position angles were fixed 
to 90°. The goal of these simulations is to characterize biases 
and errors; the uniformity of the parameter distributions adopted 
here is therefore perfectly suitable, even though real galaxy pa- 
rameters may not be so distributed. For the same reason, we do 
not take any redshift effect into account. Each simulation was 
convolved with the reconstructed PSF as explained in |5.2.2| The 
same PSF was used in both creating and analyzing the simula- 
tions, so the results will not include any error due to PSF mis- 
match. In order to simulate background noise, objects are em- 
bedded at random positions in the fields and detected with the 
same SExtractor parameters as for the real sources. 



Fig. 5. Le Phare photometric redshift distribution for the 60 
matched objects. The distribution is peaked around z ~ 0.8, 
in good agreement with the predictions of simple PLE models. 
Error bars show poissonnian errors. 

computed with the code Le Phar^ A standard method is 
implemented, including an iterative zero-point refinement com- 
bined with a template optimizat ion procedure and the appli- 
cation of a Bayesian approach jllbert et al.l l2006ah . We used 
1095 spe ctroscopic redshifts taken from the zCOSMOS Survey 
jLillv & The Zcosmos Team. 2006) to measure the photometric 
redshifts. This method allows to reach an accuracy of cr^-/{l + 
Zs) - 0.031 with 1.0% catastrophic errors, defined as Az/(1 -i- 
z.)>0.15. 

The multi-color catalog of the COSMOS survey (Capack 
et al. 2006) consists of photometry measurements over 3 arc- 
second diameter apertures for deep B j, Vj, g+, r+, i+, z+ Subaru 
data taken with SuprimeCam, m*, ;* bands with MegaCam 
(CFHT), u',g',r',i',z' information from the Sloan Digital 
Sky Survey (SDSS), K, magnitude from KPNO/CTIO, and 
F816W HST/ACS magnitude. Objects were matched between 
the COSMOS and the NACO catalogs within a radius of 2" 
which takes possible astrometry differences between the cala- 
logues into account. We matched 60 objects out of the 79 de- 
tected in the NACO fields. Fig. [5] shows the photometric red- 
shift distribution for these 60 matched objects. As expected for a 
galaxy sample l imited to Ks = 22, th e redshift distribution peaks 
around z ~ 0.8 (iMignoli et al.ll2005h . 

5. Automated morphology classifications 

The 79 objects identified as extended sources are morpholog- 
ically classified using two automated methods based on direct 
model fitting and on learning classification. Automated classifi- 
cations have the fundamental properties of being objective, thus 
reproducible, and they allow a precise error characterization. 
We proceeded in two steps; first we detected irregular objects 

^ Ihttp : //www . lam . oamp . £r/arnouts/LE_PHARE . htmll 



5.1. Irregular objects 

The detection of objects presenting irregulariti es is made us- 
ing the concentration and as ymmetry estimators dAbraham et al.l 
119941: [Abraham et al.1 [19961) . Concentration is computed as the 
ratio of the flux between the inner and outer isophotes of nor- 
malized radii 0.3 and 1 within the isophotal area enclosed by pix- 
els 3cr above the sky level. The corresponding limiting surface 
brightness varies between fi = 18.76 - 20.39 mag.arcsec^^ be- 
cause of the variations in the exposure time between the different 
fields and the intrinsic variation of the sky level, which is impor- 
tant for infrared ground-based observations. Indeed, the value of 
C is quite sensitive to the estimate of the background level since 
an error in this value will result in different limiting isophotes, 
and a fraction of the galaxy flux can be lost. To estimate the er- 
ror in C introduced by the different isophote levels, we computed 
the variations in the C value for variations in the limiting surface 
brightness of Afi - I. We found a fairly small error (AC ~ 0.06), 
so we decided not to apply any corrections. Asymmetry (A) was 
obtained by rotating the galaxy image about its center by 1 80 and 
self-substracting it to the unrotated image after sky substraction. 
Local sky level is estimated using SExtractor output parameter 
BACKGROUND that gives the background level at the galaxy 
centroid position. The center of rotation was determined by first 
smoothing the galaxy image with a Gaussian kernel of cr = 1 and 
then choosing t he location of the max imum pixel as the center, 
as explained in [Abraham et al.l (Il996h . Since the absolute value 
for the residual light is used, noise in the images shows up as a 
small positive A signa l even, in perfectly symmetrical objects. 
(IConselice et al.ll2000h . This is why we applied a noise correc- 
tion to the computation: the value of A in a portion of sky with 
area equal to that enclosed by the galaxy isophote.The definition 
of A that is used in the rest of the paper is: 

_ 1 / Z m j) - /i8o(/, ;)l Z \B(i, j) - gisoQ-, 
2 1 Zmj) 2 /(/,;■) / 

To establish the boundaries between regular and irregular ob- 
jects, a calibration of the C-A plane is needed. We thus simulated 
a set of galaxies with different galaxy parameters and magni- 
tudes ranging between 17 < Ks < 23 that we embedded in the 
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Fig. 7. Asymmetry versus magnitude: asymmetry begins to grow 
only at magnitudes greater then 22.2 which is greater than the 
limiting magnitude. 



real images. We computed the C and A parameters of these ob- 
jects and plotted the C-A plane. Since irregular objects cannot 
be simulated in a meaningful way, a kind of extrapolation was 
employed, based on two facts: a)irregular galaxies have flatter 
photometric profiles (less concentrated) and to be more asym- 
metric than regular objects; b)those objects are not present in the 
simulated sample. We thus defined the irregular zone as the up- 
per left corner of the C-A plane where no simulated objects are 
found (Fig.|6j. 

As said previously, rotational asymmetry is affected by noise 
even after correction. This means that fainter objects might ap- 
pear more asymmetric and can thus induce a bias in the num- 
ber of irregular objects at the faint end of the sample. To esti- 
mate this error, we plotted the asymmetry parameter for a sam- 
ple of 1000 simulated galaxies with magnitudes ranging between 
17 < Ks < 23 (Fig.|7]). The plot shows that asymmetry begins 
to be aff'ected by noise only at magnitudes greater than 22.2, 
which is the magnitude limit of our working sample. In sum- 
mary, we found the location of the irregular/peculiar objects by 
simulating a set of regular galaxies and defining the peculiar area 
as the area outiside. Then, we plotted the observed data on this 
plane and count the galaxies in this peculiar area. We counted 
10 observed objects in this zone, i.e. 12% of the sample. We can 
attempt to quantify the error in this classification by consider- 
ing the regular simulated objects that fall in the irregular zone. 
This gives the fraction of regular objects that are misclassified. 
We counted 27 objects out of 1000. we conclude that 12% + 2.7 
of our sample corresponds to peculiar objects, in the magnitude 
range \1 < Ks < 22. 



5.2. Regular objects: disk dominated - bulge dominated 
Separation 

5.2.1. C-A morphology 

The positions of galaxies in the C-A plane are used to separate 
the bulge and disk-dominated galaxies as follows: 

Calibration: A calibration of the C-A plane is needed before 
classification to investigate where the objects exactly fall. The 
strategy followed in this paper is two-fold: first we draw the ir- 
regular border as defined in Section 15.11 The border between 
bulge-dominated and disk-dominated objects can be deduced in 
a more automated way thanks to the analysis of simulated galax- 
ies. We took the same 1000 galaxies as above, for which the 
morphological type is known, and draw their positions in the C- 
A plane (Fig.[8]l. The border is then defined with a c lassification 
method based on support vector machines d Vapnikll 1995i) Fl. We 
decided not to use classical boundaries employed in previous 
works because: a) those boundaries were not obtained in the K- 
band and b)we were looking for an objective method that did not 
require visual inspections. Support vector machines (SVM) non- 
linearly map their n-dimensional input space into a "high dimen- 
sional feature space". In this high-dimensional feature space, a 
linear classifier is constructed. SVM have two main parameters 
that can be changed: the kernel function and the tolerance C. The 
kernel function corresponds to the expected shape of the border, 
for instance, if the objects are distributed with a gaussian dis- 
tribution, then a gaussian kernel will be used. The adjustment 
of the border will also take a tolerance factor C into account. 
If C is high, the machine will not allow any object to be on 
the wrong side of the border. As a consequence, if the objects 
are strongly mixed in a given plane, the border can have a very 
complex shape. In contrast, if C is too low the machine will not 
reach an optimal separation. In a first approach, a linear kernel 
was used, assuming that the two families of objects can be sepa- 
rated with a linear function, in order to be coherent with previous 
works. Thus the C parameter is set to be infinite. The results of 
this separation are displayed in Fig. [8] 

Accuracy: Automated classifications are useful because they al- 
low a characterization of errors. Once the boundaries are drawn, 
we generated another set of 500 fake objects with known mor- 
phology that we place again in the C-A plane and that we clas- 
sifed in the three morphological types we have defined. We then 
compared the results of our classification scheme with the initial 
morphological type. Errors were estimated in magnitude bins, 
from Ks = IS to Ks = 22 (Table We achieved 70% of good 
identifications up to Ks = 21 (AB - 23) and 66% up to our mag- 
nitude limit Ks = 22. That means that we are able to classify 
galaxies in the three main morphological types with a reliable 
accuracy at least up to Kab - 23. 

Some words about the C-A plane: At first sight, the distri- 
bution of galaxies in our C-A plane looks significantly different 
that wh at has been reported in previous works using this tech- 
niques. (lAbraham et al.lll996l:lBrinchmann et al.ll 19981) . Indeed, 
the slope of the separating border between bulge and disk- 
dominated galaxies has been found to be positive, whereas the 
one found here is negative, although previous classifications are 
somewhat arbitrary. As a consequence, bulge-dominated galax- 



* |http://www.isis.ecs.soton.ac.uk/resources/svniinfo/| 




Fig. 6. Separation between regular and irregular objects. Left: simulated objects (empty squares). Right: real objects (crosses) . 
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Fig. 8. C-A calibration and classification. Boundaries are drawn using an automated classification method (S VM) that avoids the 
use of a nearby sample and subjective visual classifications. Left: simulated objects, open squares: objects with B/T < 0.2, filled 
squares: objects with B/T > 0.8. Right: real objects. . 



Magnitude 


Correct Identifications 


K, < 19 


80% 


K, < 20 


73% 


K, < 21 


70% 


K, < 22 


66% 



Table 2. Error estimates of the C-A classification. Fraction of 
misclassified objects for several magnitude ranges. 



ies lie in the top right corner of the plane in our classification 
rather than in the bottom right corner in most other studies. There 
might be several reasons to explain this effect: 

- As NACO images are undersampled, highly concentrated 
objects only cover a few pixels, consequently, even a small 
mismatch in the determination of the rotation center can lead 
to large asymmetries. To check this effect we computed con- 



centration and asymmetry parameters for the detected point- 
like sources. We find that, indeed, they appear to have higher 
asymmetric values than the galaxies. 
- It might also be a conseq uence of the m ethod used 
for the boundary estimate. [Abraham et al.l (1 19961) and 
iBrinchmann et al.f ([l998t) used a visual inspection based on 
a local survey and, in order to account for the effects of red- 
shift, they applied corrections to the concentration parame- 
ter. In this work, we used a fully automated method based on 
simulations that reproduce the observational conditions very 
closely and on learning classification methods, so no correc- 
tion is needed in principle, but this gives a different border 
The question that consequently arises is whether the differ- 
ences in the boundaries produce important differences in the 
classification procedure. All this points are fully discussed in 
Sect. |6] thanks to detailed comparisons with space observa- 
tions of the sample. 
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Fig. 9. Example of PSF fitting in Fourier space. Squares: 
observations, dashed line: seeing-limited MTF, dotted line: 
diffraction-limited MTF. The AO MTF contains higher frequen- 
cies than the seeing-limited one. The telescope diffraction limit 
is not reached however in this example due to the undersampling 
of the instrumental setup. 



5.2.2. Model-fitting morpliology 

The second method is based on a direct two components fit- 
ting with exponenti al and de Vaucouleurs profiles, using GIM2D 
dSimard et al.ll2002l) . The 2D galaxy model used by GIM2D has 
1 1 parameters that are fitted to the real data. The most important 
ones are the total flux and the bulge fraction BIT {- for pure 
disk systems). Other parameters are the (i)disk scale radius r^, 
(ii)the disk inclination ;, (iii)the effective radius re (iv)the ellip- 
ticiy e of the bulge component, and other geometric parameters 
for the center and orientation of both components. As GIM2D 
also estimates the local sky level using image statistics, this is 
the value that is used. 



PSF reconstruction: To obtain reliable results, GIM2D needs a 
noise-free, well-sampled PSF. This is why special attention has 
been paid to PSF reconstruction. Here, classical methods, such 
as DAOPHOT or Tiny Tim, could not be used for two reasons. 
First, the Adaptive Optics PSF has a specific shape that is neither 
'seeing limited like' nor 'spatial like'. An AO system operated 
with a guide star of moderate brightness {V ~ 14 ) can only par- 
tially correct for turbulence-induced distortions. This partially- 
corrected PSF consists of two components: a difFraction-limited 
core, superimposed on a seeing halo. Second, to have a larger 
field and a better sensitivity, data were under-sampled by a fac- 
tor 2 (0.054" pixel scale, whereas 0.02" is needed to be Nyquist 
sampled). 

We developed a simple algorithm that uses field stars to gen- 
erate Nyquist-Shannon-sampled PSFs by means of a fitting pro- 
cedure in Fourier space. The process is as follows: we generate 



PSF,„(x,y) = ^ 



SRxA(jc,3;)H-(l - SR) x exp - 



where SR is the Strehl ratio, A{x,y) the bi-dimensional Airy 
function, cr the Gaussian dispersion that can be related to the 
seeing and k is a normalization factor. This artificial PSF is built 
with a Nyquist-Shannon sampling, binned by a factor 2 to reach 
the real-image pixel scale and finally Fourier-transformed to cre- 
ate a simulated MTF (power spectrum). On the other hand, for 
each observed star, its Fourier transform is fitted with the sim- 
ulated MTF. The parameters estimated that way (SR, cr, and k) 
are then used to build an estimate of the PSF with the correct 
Nyquist-Shannon sampling. 

Working in Fourier space avoids including the background 
estimate and PSF position as a fit parameter, which is particu- 
larly delicate in our case, since the FWHM is less than 2 pixels 
large. In the few cases where the fitting procedure does not con- 
verge a second Gaussian halo is added. Fig.|9]shows the result of 
the fitting for one star in the spatial frequency domain. In this pa- 
per, we do not consider variations in the PSF caused by adaptive 
optics, such as anisoplanetism, but we are working on building a 
complete model for PSF estimate for future observations. 

Error analysis: We ran GIM2D on the 79 objects with mag- 
nitudes ranging between Ks - 17 - 22, using a two compo- 
nents model and the artificial PSF built as described above. We 
used the GIM2D mode that allows use of oversampled PSFs 
to deal with undersampled data, since the PSF recovered with 
the method explained above is Nyquist sampled. The fitting con- 
verged for the whole sample, and the results are quite convincing 
in terms of residual images. (Table [3]) 

Visual inspection of the models compared to the real images 
(Table O also reveals good agreement, in particular for bright 
sources. For the faintest objects, however, it is more difficult to 
estimate the fitting accuracy. Indeed, inspection of image resid- 
uals is not a robust accuracy test, since there may still be strong 
degeneracies even when the image residuals do not show any 
features. This is why objective and systematic error characteriza- 
tion is needed. To do so we generated a sample of 1000 synthetic 
galaxies with known galaxy parameters uniformly distributed: 
< BIT < 1, < < 0.5", < < 0.5", < e < 0.7, 
0° < / < 85°. The Sersic bulge index was fixed at « - 4, and 
both bulg e and disk po s ition angles were fixed to 90°. As ex- 
plained in ISimard et al.l (|2002|) . the goal of these simulations is 
to characterize biases and error. The uniformity of the parame- 
ter distributions adopted here is therefore perfectly suitable even 
though real galaxy parameters may not be so distributed. Each 
simulation is convolved with the reconstructed PSF. The same 
PSF is used in both creating and analyzing the simulations, so the 
results will not include any error due to PSF mismatch. In order 
to simulate background noise, objects are embedded at random 
positions in the fields and detected with the same Sextractor pa- 
rameters as for the real sources. Finally, the GIM2D output files 
are processed through the same scripts to produce a catalog of 
final recovered st ructural parameters. 

Following the lSimard et al.l (l2002h procedure, we decided to 
represent errors in a set of two-dimensional maps giving system- 
atic and random errors at each position. The GIM2D parameter 
space is a complex space with 11 dimensions, so these maps 
can only offer a limited representation of the complex multidi- 
mensional error functions but makes interpretation much sim- 
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Table 4. Error analysis of the bulge fraction B/T for different recovered magnitude ranges and different bins in recovered galaxy 
size. The galaxy size is represented by the half-light radius and is distributed into 4 bins in log(r/,;). In the top left corner bright and 
small objects are found whereas faint and large objects are placed in the bottom right corner AB/T is the average difference between 
introduced and recovered values of B/T, while crAB/T is the dispersion (see text for details). N is the number of simulations used 
for each bin. 



pier The error analysis presented in this paper focuses on the 
error made on the main morphological estimator, the bulge frac- 
tion, as a function of two main parameters: apparent magni- 
tude and half light radius. Systematic errors are computed as 
the m ean dif ference between the introduced and the recovered 
value: AB/T = ^^'^^''j^ ^^^''^ and random errors as the square root 

of the variance of the difference: crAB/T - ^ 
Table m precisely shows in details the sources of errors on B/T 
as a function of galaxy magnitude and half-light radius. 

The main result after looking at the results of simulations 
is that, for objects brighter than Ks ~ 19 {AB ~ 21), the bulge 
fraction is recovered with a bias close to zero and a random error 
around 20%. This is true even for small objects (-1 < log{rhi) < 
-0.75), and it is comparable to what is obtained for the brightest 
objects in the 1-band with HST dSimard et al.ll2002l) . For fainter 
magnitudes, we can see two main effects: 

- first, an increasing bias in the bulge fraction (~ 20%), in 
particular for large objects (/og(r/,;) > -0.7 5). This is a 
well-k nown GIM2D effect for low S/N objects dSimard et alj 
I2OO2I) : the outer wings of steep surface brightness profiles, 
such as the r'^^ profile, are hidden in low S/N objects which 
artificially decreases the recovered bulge fraction. This can 
also be a consequence of errors in the sky-level estimate that 
causes that a fraction of the flux is hided by noise. 

- second, a slight increase in the random error, which becomes 
closer to 30%. This is a consequence of an increasing degen- 
eracy of the parameters space with decreasing S/N. 

But for most of the objects brighter than = 19 galaxy 
parameters can be estimated correctly (cr ~ 0.2 and b < 0.1), 
even for small objects (r/,/ < 0.3") of size comparable the limits 
of space observations (see Section|6ll. 

5.3. Results of the analysis and comparison of classifications 

We classified the galaxies into three main morphological types 
according to the fitting results. One of the main results is that 
about 12% ± 2.7% of our sample corresponds to peculiar or ir- 
regular objects (10 objects out of 79). For the rest of the sample, 
the G1M2D analysis finds 21 (~ 26%) bulge-dominated galaxies 
(B/T > 0.5) and 48 (~ 60%) disk-dominated {B/T < 0.5) while 
for the C-A classification, we find 54 (~ 67%) disk-dominated 
galaxies and 15 (~ 19%) bulge-dominated ones. 



Looking in more detail into the relaibility of the two classifi- 
cation schemes, we did a one-to-one comparison of the morpho- 
logical types assigned by the GI1VI2D analysis or the C-A one 
(Fig [Tol l: we computed the probability that a galaxy classified 
using the G1M2D classification is classified with the same mor- 
phological type by C-A. The probability was computed by divid- 
ing the number of galaxies in each morphological C-A bin by the 
total number of galaxies of the same type selected with G1M2D. 
Overall, there is good agreement between both classifications in 
the whole sample. The probability that a disk dominated galaxy 
identified by G1M2D has the same morphological type in C-A 
classification is p - 0.81, but only p - 0.30 for bulge-dominated 
galaxies including the faintest objects (Ks < 23). For irregulars, 
it is obviously p - I since the detection procedure is the same 
in both methods. 

There might be two reasons why the classifications are not 
exactly the same. First, the S/N might cause discrepancies. 
Indeed, as we show in Sect. 15.2.21 at low S/N, G1M2D tends to 
under estimate the bulge fraction. This implies that some galax- 
ies detected by G1M2D as disk-dominated are in fact detected as 
bulge-dominated by C-A. FigurefTOlshows the effect of reducing 
the limiting magnitude to Ks = 20: the fraction of objects clas- 
sified as bulge dominated by G1M2D and C-A rises up to 0.67. 
Second, it might be a problem of definition. Indeed, the mor- 
phological bins are not exactly the same in both classifications. 
In particular, objects with intermediate morphological type (i.e 
B/T ~0.5) might cause discrepancies. If we remove those ob- 
jects from the sample, 80% of the bulge-dominated objects and 
95% of the disk-dominated objects detected by GIM2D are also 
detected by C-A with the same classification. 

Either way, the comparison of both classifications allows a 
quantification of the error in classification of regular galaxies in 
the sense that it seems reasonable to think that the true value 
should be somewhere between the two results. The GIM2D es- 
timate thus gives a lower limit for the early-type fraction and C- 
A the upper-limit and vice versa for the late-type fraction. This 
way, we conclude that the mixing of population in our sample is: 
24% +4% of early-type galaxies, 64% +4% of late-type galaxies, 
and 12% + 2.7% of irregular/peculiar galaxies. 

Our results offer the first direct measurement of the distri- 
bution of galaxy in three morphological types at z ^ 1 from 
high spatial-resolution imaging in the K-band. We observe that 
the fraction of 12% of irregular objects at z = 1 is signif- 
icantly higher than the fraction of these objects in the local 
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Table 3. GIM2D output for some objects. Left: original images. 
Middle: GIM2D models. Right: Residuals. The small thumb- 
nails show the real and the model galaxy. The image size is 
1.7' X 1.7' 



GIM2D 
K<23 K<20 




Bn-^o.z B/T>o.a irr BfT<^0.2 BfT >o.B Irr 



Fig. 10. Comparison of classification methods, show the prob- 
ability that a galaxy classified with G1M2D is classified in the 
same morphological type by C-A. (see text for details). 
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Fig. 11. Example of classification in the three main morphologi- 
cal types at different redshifts. The image size is 1.7 x 1.7 



Universe, confirming from rest-frame data at ~ 1 microns the 
well docu mented trend of this pop ulation increasing with red- 
shift (e.g. iBrinchmann et al.l (Il998l) ). However, this result must 
be taken with caution. Indeed GIM2D accuracy decreases for ob- 
jects fainter than Ks = 19, which represents 80% of the sample. 
Moreover at the faint end, the fraction of irregular objects can be 
overestimated because of the low S/N. But there are good rea- 
sons to consider this result significant. Even though there is an 
over estimation of disks in the faint end, the morphological clas- 
sification bins are large enough to reduce the number of false 
classifications. Indeed, even in the zones where the random er- 
ror in the bulge fraction estimate is ~ 0.3 or larger, we do not 
classify a pure bulge (B/T ~ 1) as a disk. 



6. Comparison with ground-based and HST 
observations 

In this section we compare our AO observations with ground- 
based and space observations. 

6.1. Ground-based observations 

Effective radii of local galaxies, except for compact dwarf galax- 
ies, range from ~ 1 . to ~ 10 kpc depe nding on their luminosity 
dBender et al.lll992l : llmpey et alJl996l) . Our spatial resolution of 
~ 0.1" corresponds to about 1 kpc at z ~ 1 and we should be able 
to determine morphological types even at z > 1 . Thus, in order to 
estimate the performance of AO deep imaging and to justify the 
automated morphology classification, we compared our images 
with deep I-band seeing-limited images taken from the Canada- 
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France-Hawaii Telescope Legacy Survey (CFTHLSfl One of 
the so-called deep fields is centered on the COSMOS field, al- 
though it is smaller than the total COSMOS area (1 square de- 
gree out of 2). Here we used the release T0003 images (March 
2006jJ especially the deep /' one, corresponding to a total inte- 
gration time of 20 hours, with an average FWHM of ~ 0.7". 

We compared real data by selecting a galaxy classified as a 
disk by GIM2D and C-A in the NACO data and by comparing 
it to the results obtained with CFHTLS data. We computed the 
surface brightness profile within the isophotal area enclosed by 
pixels 3cr above the sky level. The corresponding limiting sur- 
face brightness is = 20 mag.arcsec^^ for the NACO image 
and lu - 25 mag.arcsec^^ for the MegaCam image. Sky levels 
and the corresponding isophotal areas were both determined us- 
ing SEXTRACTOR. 

The surface brightness profile was fitted with both a PSF- 
convolved de Vaucouleurs profile and a PSF-convolved expo- 
nential profile. Figure [12] shows that, with seeing-limited obser- 
vations, it is more difficult to establish a clear separation be- 
tween both profiles at small distances from the galaxy center 
(i.e. ~ 0.2 ), even if the determination of the brightness profile is 
possible at much larger distance (i.e. ~ 1 ) thanks to the depth of 
the images and the low noise level of the sky background. This 
supports the results obtained with GIM2D, which show that a 
correct estimate of the bulge fraction is possible for small ob- 
jects. Although ultra-deep sub-arcsecond imaging is powerful in 
terms of high number statistics, thanks to the wide field cover- 
age, we consider that it is more rewarding to look at a smaller 
sample of galaxies, but with more reliable morphology determi- 
nations thanks to the spatial gain of the AO. 

6.2. Space observations 

We compare our images with space data taken from the 
COSMOS survey. Since our observed fields were selected within 
the COSMOS area, the same objects were observed with the 
HST-ACS in the I-band at high spatial resolution. We thus mor- 
phologically classified the 60 objects for which the photometric 
redshift are known (Sect. HJi. We used those results to both es- 
timate the effect of the observation band on morphology and to 
validate our method to divide the C/A plane. The C-A estima- 
tors were calibrated with simulations using the same method as 
for the K band data. Standard boundaries, from other existing 
works, were used to divide the C-A plane. Figure [16] shows the 
C-A plane cut. The figure also shows the border between bulge- 
dominated and disk-dominated galaxies obtained with the au- 
tomatic method described in Sect. |5.2.Tl for this population. We 
again find a negative slope for the border between disk and bulge 
dominated objects. We find for the whole sample 32% + 1 .6% ir- 
regulars, 47% + 1.5% disk-dominated, and 21% + 2.5% bulge 
dominated. 



6.2.1. About boundaries 

As said, the computed boundaries of the C-A plane are different 
from what can be found in the literature. Previ ous works have 
been done in the 1-band usin g HST imaging ([Abraham et alj 
119961: iBrinchmann et al.|[T998l) . As we have a sample observed 
in the I-band, we are able to establish whether the change in the 
boundaries has significant consequences in the morphological 
classifications. To do so, we classified the 1-band sample using 
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Fig. 13. Galaxy distribution: comparison between K-band and I 
band C-A classifications. The figure shows the probability that 
a galaxy in the K-band is classified in the same morphological 
type in the 1-band. 
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the IBrinchmann et al.l (1 19981) boundaries and compared the re- 
sults to the ones obtained with our method (Fig. [T4] |. We find 
that there are no significant discrepancies between both classifi- 
cations. We conclude that our method is valid and moreover has 
the key advantage being free of subjective judgments. 

6.2.2. Rest-frame morphologies 

We observe some discrepancies in the global morphological dis- 
tributions between the 1 and K bands, in particular more per- 
turbed morphologies are seen in the I band. When we look 
at each object individually (Fig. [TJb we confirm this trend: 
there are uncertainties between K irregulars and / disks and 
between / disks and K bulges. Indeed an important fraction of 
bulge-dominated objects and disk-dominated objects detected by 
NACO are seen as disk-dominated and irregulars respectively, 
by ACS, as if there was a systematic trend that moves objects to 
later types when we move to shorter wavelengths. Certainly, the 
number of objects is small and a few mismatches cause high dis- 
crepancies in Fig. [13] However, this is an expected effect since 
ACS probes younger stellar populations. A visual inspection of 
some of the objects that present different morphologies reveals 
that some of the ACS irregulars are in fact well-resolved spiral 
galaxies with inhomogeneities that probably increase the asym- 
metry indices. 

In order to correctly compare both classifications we need to 
correct the measurements to estimate how galaxies would look 
if they were obse rved locally in the same p hotometric band. As 
a matter of fact, IBrinchmann et alj ([1998) showed that high-z 
galaxies imaged by HST differ in appearance from their local 
counterparts because of their reduced apparent size and sampling 
characteristics, a lower S/N and reduced surface brightness with 
respect to the sky background and a shift in the rest wavelength 
of the observations. These effects combine to give some uncer- 
tainty in the morphological classifications of galaxies. 
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Fig. 12. Comparison with ground-based observations. We performed a profile fitting on the same real galaxy observed with NACO 
(left) and MegaCam (CFHTLS-;' band, right). The galaxy magnitude is Kab = 20.5 and /' = 21.3. Surface brightness profiles were 
computed within the isophotal areas enclosed by pixels 3cr above the sky level. The corresponding limiting surface brightness is 
yu = 20 mag.arcsec'^ for the NACO image and ji - 25 mag.arcsec^^ for the MegaCam image. The fit was done with a pure de 
Vaucouleurs and exponential profile. 



The first effect is a chan ge in the concentration v alue mea- 
sured at low redshift. Indeed. lBrinchmann et alj d 19981) draw the 
boundaries in the C-A plane using a local sample ( iFrei et alj 
1 19961) visually classified. However, the concentration value de- 
pends on redshift, since the threshold is defined relative to the 
sky. Thus, less of the galaxy is sampled because of cosmological 
dimming. The solution they adopt is to correct C for this effect. 
We do not need a correction of the concentration in this paper 
because we use simulations that reproduce exactly the observing 
conditions to calibrate the C-A plane. The result is that bound- 
aries are moved with respect to a local classification instead of 
changing the C value. 

The shift in the rest-frame wavelength of observations is 
however more difficult to estimate. Indeed the question that 
arises here is whether the morphological type estimated at high 
redshift would be the same if observed at low redshift. When 
observing a galaxy in the K-band at redshift z ~ I, the equiva- 
lent rest-frame wavelength will be around the z band, however, 
when observed in the I-band, the rest-frame will be around the 
B band. That implies that a mismatch can exist in the morpho- 
logical classification since we are not probing the same mor- 
phological blocks. To correct for this effects we need to 'move 
the objects' into a common rest-frame wavelength. This is the 
called morpholog i cal k correction. The method employed by 
iBrinchmann et al.l (Il998h consists in determining the morphol- 
ogy from a local sample, redshifting the objects using SED mod- 
els, and looking at the fraction of galaxies that move in to an 
other morphological class. A drift coefficient that characterizes 
the drift from category X to category Y is thus defined as 



Once the fraction of missclassified objects is determined, the ob- 
served number of objects in class X can be related to the true 
number: 



Here we proceed as follows i fBrinchmann et al.l ( Il998l) com- 
puted the coefficients to shift from the I observed morphology to 
the R rest frame morphology, we use those coefficients to correct 
the observed HST morphology of our sample to the one observed 
in the R rest frame band, since the filter used for observations is 
the same (F814W). Once we have this corrected morphology, we 
can compare it to the NACO uncorrected morphology. This can 
be done because the observed sample is exactly the same in the 
K and in the I-band. If we were in the same rest-frame band, we 
should find the same morphology. 

We use the coefficients computed by IBrinchmann et al.l 
(1 19981) to correct the ACS morphology and divide the sample 
into two redshift bins (z < 0.8 and z > 0.8). Results are shown 
inPig.ini 



- At low redshift (z < 0.8) the corrected I-band distributions 
and the uncorrected K-band are similar The fact that there 
is no significant difference between the rest-frame R band 
and the uncorrected K-band population indicates that, when 
observing in the K-band, there is no need for morphological 
k correction. 

- At higher redshift, the K-band distribution tends to give a 
larger fraction of bulge-dominated objects and a lower frac- 
tion of irregular objects than the I-band data. This can be 
caused by the fac t that the drift coefficients computed by 
IBrinchmann et al.l (Il998l) might contain errors. Since they 
are based on SED fitting, it is logical to think that errors are 
more important for high-redshift objects. This would explain 
why the results are consistent at low redshift. It can also be 
a consequence of the size of the sample: as the number of 
objects per bin is small, errors in the classification of a small 
fraction can lead to discrepancies between the two bands. 
This issue will certainly be solved with a larger sample. 
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Fig. 16. C-A cut for the ACS images. The same classification procedure has been repeated for the same sample observed with ACS 
in the I-band. Left: simulated o bjects. Right: real object s. Circles: BjT < 0.2, Filled squares: BjT > 0.8, Crosses: real objects. 
Dotted line is the border used in lBrinchmann et aTl (Il998l) to separate bulge dominated from disk dominated, dashed Une is the one 
computed in the paper. 
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Fig. 14. Comparison of classifications with different boundaries. 
We repeat the morphological classification with the boundaries 
used by [Brinchma nn et al. ( 1998). We conclude that the results 
do not change significantly which supports the validity of the 
employed method. 



7. Summary and conclusions 

We analyzed the morphologies of a sample of 79 galaxies in the 
near-infrared thanks to adaptive optics imaging at a resolution of 
0.1 . Thanks to extensive simulations, we showed that adaptive 
optics can be used to obtain reliable high-resolution morphologi- 
cal information in an automated way and is thus adapted to large 
observation programmes: 



- Galaxy parameters (bulge fraction, disk scale length, bulge 
effective radius) can be estimated by means of model fitting 
with errors lower than 20% up to Ks - 19. This is compara- 
ble with space data. 

- Galaxies can be separated into three main morphological 
bins up to Ks -21, with at least 70% of good identifications 
and 66% up to K,, = 22. 

- For fainter galaxies (22 < Ks < 23) the accuracy of the mor- 
phological classification decreases. Simulations show that 
bulge-dominated and disk-dominated galaxies can be sepa- 
rated with only 55% accuracy at Kg = 23. 

We obtain, for the first time, an estimate of the mix of 
morphological types of the galaxy population up to z ^ 1 
from ground-based K-band observations with high spatial res- 
olution comparable or better than visible imaging from space. 
We demonstrate that estimating morphology from K-band data 
at z ^ 1 is not affected by morphological k correction, as there 
is no significant difference between our population and the cor- 
rected I-band population. We find that the fraction of irregulars 
at z ^ 1 is about 12% + 2.7% using automated classification 
methods. This is higher than what is found in local surveys, con- 
firming the well-established trend toward an increasing fraction 
of irregular galaxies with redshift as observed from surveys in 
the visible. Our small sample does not allow us to reach firm 
conclusions on the evolution of the fraction of late-type or early- 
type galaxies, but classifying galaxies from K-band AO imaging 
data is demonstrated to be reliable. 

From this work it seems clear that adaptive optics can be 
used for observational cosmology with reliable accuracy, and 
consequently data of this type should contribute to a better un- 
derstanding of galaxy evolution in the future. However, it is still 
a new technique and technical difficulties exist, such as variable 
PSF, small fields, subsampling and the need of guide stars that 
make the use of classical reduction methods more difficult. This 
is now easier with laser guide stars becoming available and new 
sets of utilities like the ones we are developing to enable easy 
data processing and analysis of adaptive optics data for the com- 
munity. This opens the way to observing the large samples re- 
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Fig. 15. Redshi ft distribution for the three morphological types. 
We plotted the iBrinchmann et al.l (1 19981) sample (circles) and 
our sa mple observed with ACS ( squares) and with NACO (tri- 
angles). |BrinchmannetaL| (Il998h and ACS data are corrected to 
the R rest-frame band. The NACO sample is observed from the 
K-band and no correction has been applied. The ACS and NACO 
samples have been separated into two redshift bins {z < 0.8 and 
z > 0.8). The represented redshifts are the median redshifts of 
each bin. 



quired to reach a robust statistical accuracy. We are planning to 
enlarge our sample by observing a large number of areas around 
bright stars in the COSMOS field, which will strongly reduce 
uncertainties in the study of morphological evolution. 

Acknowledgements. The authors want to thank the anonymous referee for many 
useful comments that greatly improved this paper 
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Table 6. Summary of the morphological classifications for the 79 detected objects. For each object we show the I and K band magnitudes and 
the estimated morphological type from AO imaging in the K-band (with GIM2D and C-A) and from HST-ACS in the I-band. BD stands for 
bulge-dominated, DD for disk-dominated and I for irregular. 



Object ID 

NHzG J100016.4+021555 
NHzGJ100016.3+021643 
NHzGJ100017.2+021637 
NHzGJ100016.3+021631 
NHzG J100015.0+021629 
NHzGJ100014.7+021630 
NHzGJ100014.8+021629 
NHzGJ100014.8+021629 
NHzG J100015.9+021629 
NHzG J100017.2+021624 
NHzG J100017. 3+021620 
NHzGJ100014.8+021624 
NHzGJ100016.7+021618 
NHzG J100016.7+021610 
NHzG J100014.9+021607 
NHzGJ100014.9+021603 
NHzG J100017. 3+021601 
NHzGJ100053.2+021934 
NHzG J100051.9+022011 
NHzG J100052.0+022008 
NHzG J 10005 1.8 +022006 
NHzG J 100050. 8 +022002 
NHzG J 100050. 8 +022000 
NHzGJ100051.9+021944 
NHzGJ100051.0+021942 
NHzGJ100052.4+021941 
NHzG J 1000 11. 1+020629 
NHzG JlOOOlO.0+020624 
NHzG JlOOOl 1.6+020617 
NHzGJlOOOlO.2+020612 
NHzGJlOOOlO.6+020612 
NHzG J100009.7+020610 
NHzG J09595 1.4+0205 14 
NHzG J095953.9+020507 
NHzG J095952.3+020459 
NHzG J095952.3+020458 
NHzG J095952.3+020448 
NHzG J095953.0+020446 
NHzG J09595 1.5+020444 
NHzG J100013. 8+020838 
NHzG J 1000 15. 3 +020923 
NHzG J100013.9+020920 
NHzGJ100014.0+020919 
NHzGJ100014.1+020918 
NHzG J100015. 1+020915 
NHzG J100016.3+020912 
NHzGJ100013.9+020909 
NHzGJ100013.9+020903 
NHzG J 1000 15. 4+020854 
NHzG J100016.1+020854 
NHzG J 100003. 4+0207 11 
NHzG J100003.6+020706 
NHzG JlOOOOl.7+020704 
NHzG J100002.5+020701 
NHzG JlOOOOO.9+020701 
NHzG J100003.2+020700 
NHzG J100002.6+020659 
NHzG JlOOOOO.9+020652 
NHzG J100002.1+020650 



RA 


DEC 


Ks 


I 


rh(arcsec) 


ZPHOT 


G2D(K) 


C-A(K) 


C-A(I) 


150.069 


2.26541 


20.1822 


99.5000 


0.296730 


99.9000 


DD 


DD 


N/A 


150.068 


2.27872 


20.7045 


22.4291 


0.333504 


0.711700 


DD 


DD 


I 


150.072 


2.27703 


22.8723 


24.1406 


0.205416 


0.247200 


DD 


DD 


DD 


150.068 


2.27552 


19.7637 


24.2207 


0.566514 


1.78470 


DD 


DD 


I 


150.063 


2.27500 


18.7053 


99.5000 


0.489564 


99.9000 


DD 


DD 


N/A 


150.061 


2.27509 


22.1654 


23.3978 


0.209682 


0.832200 


BD 


DD 


DD 


150.062 


2.27487 


22.0163 


23.7523 


0.169614 


0.793700 


I 


I 


DD 


150.062 


2.27478 


21.3202 


99.5000 


0.292356 


99.9000 


I 


I 


N/A 


150.067 


2.27476 


21.7583 


22.6455 


0.246294 


0.0400000 


BD 


DD 


DD 


150.072 


2.27355 


20.0909 


23.4449 


0.243324 


0.956100 


BD 


DD 


DD 


150.072 


2.27223 


20.2320 


99.5000 


0.403920 


99.9000 


I 


I 


N/A 


150.062 


2.27360 


19.8198 


99.5000 


0.231174 


99.9000 


BD 


BD 


N/A 


150.070 


2.27174 


19.0273 


20.9238 


0.640170 


0.219100 


DD 


DD 


DD 


150.070 


2.26963 


20.3879 


22.7999 


0.271026 


0.740000 


BD 


DD 


BD 


150.062 


2.26870 


18.5753 


21.8364 


0.392418 


0.875600 


BD 


BD 


DD 


150.062 


2.26755 


22.2439 


22.7490 


0.235656 


0.475700 


BD 


DD 


BD 


150.072 


2.26712 


20.0121 


23.8199 


0.266760 


1.33790 


DD 


DD 


I 


150.222 


2.32621 


18.5348 


21.3553 


0.531684 


0.886100 


DD 


DD 


DD 


150.216 


2.33658 


22.6880 


22.6407 


0.183600 


0.810100 


BD 


DD 


DD 


150.217 


2.33579 


19.2141 


21.9404 


0.250614 


0.832600 


DD 


BD 


BD 


150.216 


2.33511 


21.9112 


23.6353 


0.254880 


0.818800 


DD 


DD 


I 


150.212 


2.33398 


20.2005 


21.8439 


0.471420 


0.514000 


DD 


DD 


DD 


150.212 


2.33337 


17.6825 


19.9774 


0.690012 


0.400600 


DD 


DD 


DD 


150.217 


2.32890 


18.7411 


20.8672 


0.396252 


0.670000 


DD 


DD 


DD 


150.213 


2.32853 


18.3509 


20.8194 


0.614520 


0.739400 


DD 


DD 


DD 


150.219 


2.32814 


20.8822 


22.7814 


0.404622 


0.799500 


I 


I 


DD 


150.047 


2.10808 


17.7249 


19.9124 


0.348300 


0.506000 


BD 


BD 


DD 


150.042 


2.10673 


19.5877 


22.7680 


0.200826 


0.987200 


DD 


BD 


DD 


150.048 


2.10486 


19.5207 


21.0617 


0.346680 


0.270200 


DD 


DD 


DD 


150.043 


2.10360 


20.3156 


99.5000 


0.155034 


99.9000 


DD 


DD 


N/A 


150.045 


2.10348 


20.5515 


99.5000 


0.167130 


99.9000 


BD 


DD 


N/A 


150.041 


2.10279 


21.2964 


99.5000 


0.178092 


99.9000 


DD 


DD 


N/A 


149.964 


2.08732 


20.8670 


21.1830 


0.379188 


0.671400 


I 


I 


DD 


149.975 


2.08549 


21.8739 


23.2931 


0.170046 


0.703700 


DD 


DD 


DD 


149.968 


2.08319 


19.9646 


99.5000 


0.358668 


99.9000 


DD 


DD 


N/A 


149.968 


2.08286 


21.7629 


99.5000 


0.202284 


99.9000 


BD 


DD 


N/A 


149.968 


2.08006 


21.3778 


99.5000 


0.178254 


99.9000 


BD 


DD 


N/A 


149.971 


2.07970 


20.2794 


23.9269 


0.192024 


0.754600 


DD 


DD 


I 


149.965 


2.07890 


21.1099 


23.6728 


0.160164 


1.37440 


DD 


BD 


DD 


150.058 


2.14397 


22.4294 


22.7827 


0.185220 


0.240300 


DD 


BD 


DD 


150.064 


2.15644 


19.5222 


21.1499 


0.395604 


0.979400 


BD 


BD 


DD 


150.058 


2.15569 


21.0909 


99.5000 


0.189000 


99.9000 


BD 


DD 


N/A 


150.058 


2.15551 


21.2097 


23.5849 


0.247158 


0.994000 


BD 


DD 


DD 


150.059 


2.15517 


19.5826 


23.5148 


0.193482 


1.26360 


BD 


BD 


DD 


150.063 


2.15427 


23.0568 


99.5000 


0.139968 


99.9000 


BD 


DD 


N/A 


150.068 


2.15348 


21.6753 


23.0442 


0.156600 


0.935900 


DD 


DD 


DD 


150.058 


2.15258 


21.4149 


22.7981 


0.306612 


1.08450 


DD 


DD 


DD 


150.058 


2.15105 


21.4370 


23.5968 


0.185112 


1.11800 


BD 


DD 


DD 


150.064 


2.14850 


20.9069 


23.9274 


0.283014 


1.16300 


DD 


DD 


I 


150.067 


2.14840 


22.1741 


99.5000 


0.168156 


99.9000 


BD 


DD 


N/A 


150.015 


2.11984 


20.5216 


22.3691 


0.432864 


0.714900 


DD 


DD 


DD 


150.015 


2.11842 


19.0519 


22.4869 


0.408024 


1.17460 


DD 


BD 


DD 


150.007 


2.11780 


17.0060 


18.8229 


0.862596 


0.338100 


DD 


DD 


DD 


150.011 


2.11709 


20.4333 


99.5000 


0.220644 


99.9000 


DD 


DD 


N/A 


150.004 


2.11697 


21.7613 


22.4334 


0.235710 


0.787500 


I 


I 


DD 


150.013 


2.11682 


18.8483 


22.2691 


0.440532 


1.06140 


DD 


DD 


I 


150.011 


2.11655 


17.6231 


99.5000 


0.397062 


99.9000 


DD 


BD 


N/A 


150.004 


2.11461 


19.4302 


20.6107 


0.372168 


0.698200 


BD 


DD 


DD 


150.009 


2.11392 


18.8301 


21.3752 


0.445932 


0.689900 


DD 


DD 


DD 
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Table 6. continued. 



Object ID 


RA 


DEC 


Ks 


I 


rh(arcsec) 


ZPHOT 


G2D(K) 


C-A(K) 


C-A(I) 


NHzG J 1 00003 . 1 +020648 


150.013 


2.11341 


20.3555 


21.9584 


0.402084 


0.422300 


DD 


DD 


I 


NHzG J 1 00003 . 1 +020642 


150.013 


2.11180 


21.3786 


23.4545 


0.283662 


0.770500 


I 


I 


DD 


NHzG J 100002. 2 +020641 


150.009 


2.11151 


20.8944 


21.5886 


0.404406 


0.334700 


I 


I 


DD 


NHzG J100002. 1+020635 


150.009 


2.10981 


21.4976 


99.5000 


0.294678 


99.9000 


DD 


DD 


N/A 


NHzG J095956. 8+020423 


149.987 


2.07308 


19.8462 


23.6537 


0.332424 


1.00700 


DD 


DD 


DD 


NHzG J095956. 7+020421 


149.987 


2.07260 


22.6499 


22.7361 


0.321570 


1.10000 


DD 


DD 


I 


NHzG J095955. 6+020420 


149.982 


2.07236 


20.8560 


22.4354 


0.202986 


0.689900 


DD 


BD 


BD 


NHzG J095955. 6+020419 


149.982 


2.07204 


21.3234 


22.7710 


0.204984 


0.624900 


DD 


DD 


BD 


NHzG J095956.4+020416 


149.985 


2.07137 


20.8485 


22.3323 


0.404082 


0.0400000 


DD 


DD 


BD 


NHzG J095956.4+020412 


149.985 


2.07027 


19.3325 


99.5000 


0.389070 


99.9000 


DD 


BD 


N/A 


NHzG J095955.8+020412 


149.983 


2.07019 


22.9091 


23.4010 


0.194778 


0.525700 


I 


I 


BD 


NHzG J095956. 8+020406 


149.987 


2.06838 


20.3433 


23.9523 


0.259470 


1.54810 


BD 


BD 


DD 


NHzG J095955. 8 +020407 


149.983 


2.06874 


18.9586 


22.0782 


0.264924 


1.20850 


DD 


BD 


DD 


NHzG J095955. 7+020404 


149.982 


2.06801 


20.4066 


23.8483 


0.391338 


1.48830 


I 


I 


DD 


NHzG J095956.2+020401 


149.984 


2.06707 


21.2254 


99.5000 


0.249372 


99.9000 


DD 


DD 


N/A 


NHzG J095955. 8+020358 


149.983 


2.06630 


21.2308 


23.8229 


0.434916 


0.881200 


DD 


DD 


I 


NHzG J095956.5+020355 


149.986 


2.06550 


18.0002 


21.4034 


0.622458 


0.794100 


DD 


DD 


DD 


NHzG J095956.2+020353 


149.984 


2.06479 


19.3081 


21.7327 


0.394146 


0.668900 


DD 


DD 


DD 


NHzG J095956.4+020353 


149.985 


2.06474 


22.1424 


23.8900 


0.188568 


0.743900 


DD 


DD 


DD 


NHzG J095957.2+020347 


149.989 


2.06312 


21.8247 


23.8017 


0.184248 


0.214600 


DD 


DD 


DD 



18 



Huertas-Company et al: Morphology of z ~ 1 galaxies with adaptive optics 



List of Objects 



'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 
'NHzG 



J100016 
J100016 
J100017, 
J100016, 
J100015, 
J100014, 
J100014, 
J100014, 
J100015, 
J100017, 
J100017, 
J100014, 
J100016, 
J100016, 
J100014, 
J100014, 
J100017, 
J100053, 
J100051, 
J100052, 
J100051, 
J 100050, 
J 100050, 
J100051, 
J100051, 
J100052, 
JlOOOll, 
JlOOOlO, 
JlOOOll, 
JlOOOlO, 
JlOOOlO, 
J100009, 
J095951, 
J095953, 
J095952, 
J095952, 
J095952, 
J095953, 
J095951, 
J100013, 
J100015, 
J100013, 
J100014, 
J100014, 
J100015, 
J100016, 
J100013, 
J100013, 
J100015, 
J100016, 
J 100003, 
J 100003, 
JlOOOOl, 
J100002, 
J 100000, 
J 100003, 
J100002, 
J 100000, 
J100002, 
J 100003, 
J 100003, 
J100002, 



4+021555 
,3+021643 
2+021637 
3+021631 
0+021629 
7+021630 
8+021629 
8+021629 
9+021629 
2+021624 
3+021620 
8+021624 
7+021618 
7+021610 
9+021607 
9+021603 
3+021601 
2+021934 
9+022011 
0+022008 
8+022006 
8+022002 
8+022000 
9+021944 
0+021942 
4+021941 
1+020629 
0+020624 
6+020617 
2+020612 
6+020612 
7+020610 
4+020514 
9+020507 
3+020459 
3+020458 
3+020448 
0+020446 
5+020444 
8+020838 
3+020923 
9+020920 
0+020919 
1+020918 
1+020915 
3+020912 
9+020909 
9+020903 
4+020854 
1+020854 
4+020711 
6+020706 
7+020704 
5+020701 
9+020701 
2+020700 
6+020659 
9+020652 
1+020650 
1+020648 
1+020642 
2+020641 



on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 
on page 



M 

m 
m 
m 



m 
m 
m 
m 
m 



m 
m 
m 
m 
m 



m 
m 
m 
m 
m 



m 
m 
m 
m 



m 
m 
m 
m 
m 



m 
m 



Huertas-Company et al: Morphology of z ~ 1 galaxies with adaptive optics 



'NHzG 


J 100002 


1+020635' 


on 


paj 


;e 


17 


'NHzG 


J095956 


8+020423' 


on 


paj 


;e 


17 


'NHzG 


J095956 


7+020421' 


on 


pa| 


;e 


17 


'NHzG 


J095955 


6+020420' 


on 


paj 


;e 


17 


'NHzG 


J095955 


6+020419' 


on 


paj 


;e 


17 


'NHzG 


J095956 


4+020416' 


on 


paj 


;e 


17 


'NHzG 


J095956 


4+020412' 


on 


paj 


;e 


17 


'NHzG 


J095955 


8+020412' 


on 


pa| 


;e 


17 


'NHzG 


J095956 


8+020406' 


on 


pa| 


;e 


17 


'NHzG 


J095955 


8+020407' 


on 


m 


;e 


17 


'NHzG 


J095955 


7+020404' 


on 




;e 


17 


'NHzG 


J095956 


2+020401' 


on 


paj 


;e 


17 


'NHzG 


J095955 


8+020358' 


on 


pa| 


;e 


17 


'NHzG 


J095956 


5+020355' 


on 


m 


;e 


17 


'NHzG 


J095956 


2+020353' 


on 


m 


;e 


17 


'NHzG 


J095956 


4+020353' 


on 


m 


;e 


17 


'NHzG 


J095957 


2+020347' 


on 


m 


;e 


17 



